Observing mergers of non-spinning black-hole binaries 
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Advances in the field of numerical relativity now make it possible to calculate t he final, most 
powerful merger phase of binary black-hole coalescence for generic binaries. The state of the art has 
advanced well beyond the equal-mass case into the unequal-mass and spinning regions of parameter 
space. We present a study of the nonspinning portion of parameter space, primarily using an 
analytic waveform model timed to available numerical data, with an emphasis on observational 
implications. We investigate the impact of varied mass ratio on merger signal-to-noise ratios (SNRs) 
for several detectors, and compare our results with expectations from the test-mass limit. We note a 
striking similarity of the waveform phasing of the merger waveform across the available mass ratios. 
Motivated by this, we calculate the match between our equal-mass and 4:1 mass-ratio waveforms 
during the merger as a function of location on the source sky, using a new formalism for the match 
that accounts for higher harmonics. This is an indicator of the amount of degeneracy in mass ratio 
for mergers of moderate mass ratio systems. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, 04.80.Nn 95.30.Sf, 95.55.Ym 97,60.Lf 


I. INTRODUCTION 

The merger of a black-hole binary will be a more lumi- 
nous source of gravitational radiation than the total elec- 
tromagnetic radiation from all the stars in the visible uni- 
verse, and will be one of the strongest sources of gravita- 
tional waves. Ground-based detectors like LIGO, Virgo, 
and GEO, currently entering their second generation of 
development, are sensitive to the mergers of stellar black 
holes, while the space-based LISA will observe mergers of 
massive and supermassive black holes. It has long been 
expected that these final mergers of black-hole binaries 
would be significant for interpreting gravitational wave 
measurements. While a physically motivated model, and 
corresponding template bank, would not be necessary for 
detection [1], such a model would be the only avenue to- 
ward extracting all of the available information about the 
system that is contained in the merger signal. Given that 
the merger is likely to constitute the majority of the de- 
tectable signal for the next generation of ground-based 
detectors 12] , such a physically motivated model would 
be necessary for gaining an understanding of the physi- 
cal sources generating the detected signals. 

In the absence of merger models, early investigations 
had to use information from perturbative approximations 
to guess at the impact of mergers. In [1], the New- 
tonian approximation for the gradual adiabatic inspiral 
of the holes, combined with the understanding of the 
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post-merger ringdown as the quasi-normal modes of a 
Kerr black hole, was used to guess at the contribution of 
mergers to the signal detectability. This guess was essen- 
tially validated by the observed behavior of numerically 
simulated merger signals [2j. However, while the power 
spectrum could be approximated, the physics behind the 
power spectrum, the amplitude and phase evolution that 
would lead to that spectrum, the accuracy with which 
the merger phase could be simulated or modeled, and 
the amount of information about the source that could 
be extracted from detected signals were completely open 
questions. Over the course of the last few years, the field 
of numerical relativity has provided a means of studying 
the detailed structure of these merger signals for the first 
time. Initially focusing only on the equal-mass, nonspin- 
ning case, several groups have since explored both the 
nonspinning axis of parameter space as well as the vast 
expanse of spinning parameter space [2-9]. The current 
availability of merger waveforms now makes it possible to 
address the questions previously mentioned through the 
measurement of these signals. 

Because the merger is the dominant contributor to the 
overall signal power, particularly for ground-based de- 
tectors where it provides the majority of the detectable 
signal, answering these questions is a critical exercise ac- 
tively being addressed by many groups. Significant atten- 
tion has been given to the problem of modeling the sig- 
nals with sufficient accuracy for detection with ground- 
based observations. Recent work has begun addressing 
not only detection, but also estimating the source param- 
eters using ground-based 110] and space-based detectors 
[11-14]. Much work has also gone into interpreting the 
available merger waveforms, in an attempt to better un- 
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derstand them, both with regard to a physical interpreta- 
tion of the source [151 and with regard to understanding 
what drives the recoil of systems due to asymmetric ra- 
diation [16— 18j. 

In this paper, we revisit the nonspinning subset of pa- 
rameter space, with the goal of studying the observa- 
tional implications of nonspinning merger waveforms. In 
Sec. II, we briefly describe the procedure for generating 

complete nonspinning waveforms. In Sec. Ill, we study 
the contribution to the achievable SNR from the inclu- 
sion of mergers, and its variation with mass ratio. In 
Sec. IV, we study more detailed comparisons of the non- 
spinning waveforms, including a novel implementation of 
the “match” statistic [19] . In Sec. V, we discuss the 
general observational implications of nonspinning merger 
waveforms. In the Appendix, we derive the formalism for 
the novel match implementation employed in Sec. IV. 


II. METHODOLOGY 

The observable quantity being measured by gravita- 
tional wave interferometers, be they ground- or space- 
based, is the strain on the spacetime, h = SL/L, or its 
derivatives. We therefore require a model of the wave- 
form, hit), that we expect to measure from black-hole 
binaries. The models we employ are predictions for the 
emitted strain, or the strain in the source frame. The 
detected strain depends on the distance to the source, 
the position on the detector's sky, and the detectors re- 
sponse. 

In this work, we focus on Advanced LIGO for ground- 
based observation, and LISA for space-based observation. 
For Advanced LIGO, we assume a constant response as 
a function of frequency, which should be adequate for all 
but the lowest-mass cases. For the detector noise, we use 
the wide-band tuning i20j typically associated with burst 
sources as was done in [2], due to the superior sensitiv- 
ity at higher frequencies where the merger will occur for 
lower masses. For LISA, we employ the effective noise 
floor from [21, 22], which includes both the contributions 
from noise sources as well as the average response of the 
detector to signals and instrumental noise, which is non- 
trivial for the higher frequencies in LISA’s band. For 
frequencies in the range 3 x 10~°Hz </< 1 x 10 -4 Hz, we 
employ a more conservative estimate of the acceleration 
noise as was done in [2], instead assuming a steeper am- 
plitude spectral density that falls off as f~ 3 [23)» Below 
3 x 1(T “Hz. we assume the detector has no sensitivity. 
We apply an overall factor of 3/20 to the LISA power 
spectral density as discussed in [24]. 

The remaining element is the model of the emitted 
waveform, which will depend on the intrinsic parameters 
(i. e. the mass, mass ratio, spin vectors, and eccentricity), 
and which will vary over the sky of the source. The 
emitted waveform can be conveniently represented by a 
harmonic mode decomposition. If h is the complex strain. 


then the mode decomposition is given by 
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where ~ 2 Ye, n are the spin- weight- (-2) spherical harmonics 
[25]. Being complex, h contains both wave polarizations, 
defined by the relationship h = h+ + ih x . For an equal- 
mass system, hit, R. 0 , o) is dominated by the quadrupole, 
the combination of 4 == 2 , m — ±2 modes: 


hquad = h.22 ^¥ 22 ( 6 , 4) + h, 2~2 2 ^2-2 (#,<£)• (2) 


Additionally, symmetry considerations for equal-mass 
nonspinning systems demand that h , 22 = h)^'- therefore 
we will often use h .22 as a proxy for the full quadrupole 
waveform. 

For the waveform comparisons presented in this work, 
we use the model first presented in [15], wdiich lias been 
validated by comparison with available data from nu- 
merical simulations for all harmonic components through 
t = 4. The model, referred to as the IRS-EOB model, 
uses the effective-one-body (EOB) Hamiltonian formal- 
ism for the inspiral [26]. For the merger-ringdown, we 
employ a novel paradigm which we call the implicit rotat- 
ing source (IRS), wherein we apply a, fit. to a physically- 
motivated functional form for the phasing (Eq. 9 in [15]), 
and. for the amplitude, a model for the flux constrained 
to be consistent, with the inspiral flux through 3.5 post- 
Newtonian (PN) order and to vanish as it approaches the 
ringdown frequency (referred to as “Model 2” and given 
by Eq. 19 in [15]). In Fig. 1, we compare the intrinsic 
error in phase and amplitude for the model and for our 
numerical waveforms, using 4:1 as a representative case. 
For brevity, we will refer to the unequal-mass runs as ra- 
tios, i, e. the q = M 1 /M 2 = 1/4 run will be the 4:1 run, 
and the ratio notation will only be used in this context. 
We shift all waveforms to peak in amplitude at t — 0, and 
to agree in phase at t = —500 M. We use the difference 
between our two highest resolutions as an indicator of our 
numerical error, and we assume errors in the model due to 
different model parameters are independent, and calcu- 
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, where A, ; is sim- 


ply A for the EOB inspiral, and A,; = {«, b. t„. Q„, Q /} 
for the IRS merger, using the notation in [15]. For S X,, 
we use the values in Table II of d5j. 

Figure 2 shows the quadrupole radiation for four mass 
ratios generated using the IRS-EOB model. The mass 
ratios shown are 1:1, 4:1, 6:1. and 20:1. We note that 
in using the 20:1 mass ratio, we have extrapolated to 
mass ratios that cannot, as yet, be validated bv simu- 
lation. The amplitudes for all the runs in Fig. 2 have 
been rescaled to better agree with the equal-mass am- 
plitude, using the leading-order Newtonian scaling. This 
also emphasizes the phasing agreement that begins in the 
late inspiral and continues through the merger waveform, 
which was discussed in J 15] and will be the topic of fur- 
ther discussion in a later section. 
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FIG. 1: Comparison of amplitude and phase errors between 
the numerical data and the IRS-EOB model, for the case of a 
4:1 waveform. The inherent phase inaccuracy of the model is 
significantly smaller than the numerical phase error, but the 
amplitude errors are comparable for the model and numerical 
data. 


Since many investigations relating to both LIGO and 
LISA have focused on detectability, rather than charac- 
terizing the signal, a model of a quadrupole-only signal 
has been adequate within their margin of error. However, 
higher harmonics can be more significant for calculations 
such as determining template fidelity with the match, or 
any attempt to extract source parameters from the sig- 
nal, as such investigations depend sensitively on the fine 
detail of the phase evolution. We will investigate the im- 
pact of higher harmonics in this context, and we will also 
include higher harmonics in our calculations of SNR, al- 
though the SNR contribution is essentially negligible for 
all cases investigated here with the possible exception of 
20:1. 


III. SNR AND POWER SCALING 

SNR is the most useful statistic for assessing the de- 
tectability of a given signal with a particular detector. 
The SNR. which we denote as p, is given by 

P 2 = (h\h), (3) 

where “(•)•)” denotes a noise- weighted inner product, 
given by 


(hilhi) = 2 df 


h*,h 2 + h x h* 2 


(4) 


and S„ (/) is the effective power spectral density of the 
detector noise 121, 22]. The sky-averaged SNR is given 
by 

<5) 

o 

Here, fi c har(/) = 2 f\h opt (f)\ is the characteristic sig- 
nal strain, and h n (f) = V5h Tms (f) = y/5fS n (f ) is the 
root-mean-square of the detector noise fluctuations mul- 
tiplied by Vo for sky-averaging. h opt (f ) is the Fourier 
transform of the optimally-oriented signal strain [1], For 
quadrupole-only cases, the sky-averaged signal strain is 
trivially calculated from the optimally-oriented strain by 
dividing by \/5. 

Before the advent of merger waveforms from numerical 
relativity, expectations about the power scaling of the 
merger waveforms, and thus the achievable SNRs, were 
formed by using the test-mass limit as a surrogate, while 
the scaling of the inspiral power can be approximated 
by PN expansions in the weak-field limit. Specifically, 
we know that the SNR from the inspiral scales as Vp to 
leading' order [ 1 ], where rj = M 1 M 2 /M 2 is the symmetric 
m.ass ratio of the binary. It was further assumed in [1], 
based on the prediction for total radiated energy in the 
test-mass limit [27], that the merger SNR scales as 77 . We 
note, however, that in [15], the peak (2,2)-mode energy 
flux was best fit by the function 

E 22 = 4.40 x 10~‘V + 5-43 x 10~ 2 r/ 4 , (6) 

with E 2 2 oc 1 / 1 , 22! 2 p 2 ■ The significant improvement 
in performance of the rj 2 + rf fit, compared to a strictly 
if fit, may indicate that differences between the physics 
of the merger for comparable masses and the test-mass 
plunge are being measured. The absence of a well-defined 
innermost stable circular orbit (ISCO) in the equal-mass 
case [28, 29], compared to the obvious ISCO for suffi- 
ciently small mass ratios, further supports this picture. 
Certainly as the masses differentiate more, the test-mass 
analogy bears out more. For the equal-mass case, Eq. 6 
indicates that the rf and r / 4 terms contribute roughly 
equally, but the rf term obviously becomes less impor- 
tant for ever-snialler mass ratios. 

The different scalings with r; are illustrated in Fig. 3, 
where we plot the Fourier transform of hybrid waveforms, 
constructed analogously to [2] by tying a PN inspiral to 
our numerical data at a point where they reach equal ac- 
curacy. We do this in part to emphasize that the change 
in scaling of the merger signal is not an artifact of the 
IRS-EOB model, but is apparent in our raw numerical 
waveforms. The left panel demonstrates the scal- 
ing of the inspiral (Mu < 0.08) for the equal-mass, 2:1, 
4:1, and 6:1 quadrupolar waveforms, and the right panel 
shows that the merger scaling is well approximated by a 
linear dependence on 77 for the merger (Mu > 0.08). The 
deviation of the peak E 22 from a simple rf scaling ap- 
pears to be due primarily to differences in the frequency 
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FIG. 2: In the top panel, we show quadrupole waveforms generated by using the model presented in [15] for mass ratios of 1:1, 
4:1, 6:1, and 20:1. When the waveforms are aligned in time based on their peak amplitudes, and aligned in phase to agree at 
said time, there is significant overlap of the waveforms for the 1:1. 4:1, and 6:1 cases over the final ~ 5 cycles leading up to 
merger (shown more clearly in the bottom panel). 


of the peak, since the signals closely following an rj scal- 
ing when evaluated at the same frequency prior to the 
peak. 

For equal masses and moderate mass ratios. 
h(t, R, 0, <j>) « hquad(f, R), so that averaging the SNR 
over the binary’s orientation is trivial. This may be suffi- 
cient for ground-based detectors if they are primarily de- 
tecting stellar-mass black holes (M < 100 Af©) due to the 
limited available mass range. If intermediate- mass black 
holes (100M© < M < 10 4 AC©) exist, then smaller mass ra- 
tios (q -C 1) may occur, and higher-order modes contain 
progressively more power relative to the quadrupole as 
q — * 0. So rather than ignoring higher-order modes or 
introducing complexity by averaging over them, we in- 
stead focus on the optimal orientation of the binary, and 
only average over the sky of the detector. 

To show the impact on SNR of varying the mass-ratio, 
we show the 1:1 and 6:1 cases for Advanced LIGO and 
the 1:1 and 20:1 cases for LISA in Figures 4 and 5, which 
demonstrate the decrease in SNR with such a significant 
deviation from equal mass. The panels show contour lines 
for both mass ratios, with one set of lines correspond- 
ing to the SNR accumulated before the corresponding' 
Schwarzsehild ISCO, and the other corresponding to the 
full signal. As described earlier, the SNR decreases as t] 
deviates from 0.25, with the inspiral SNR scaling as y/J) 
and the merger SNR scaling as r]. We therefore expect the 
SNRs to scale roughly as y/rj for lowest masses where the 
inspiral matters most, and as i] (potentially with terms 
of higher power in t) as well) for higher masses where the 


merger contributes the majority of SNR. 

In both mass-ratio cases for each detector, the late 
inspiral-merger phase constitutes the majority of the 
SNR in high-SNR events. The merger contributes sig- 
nificantly for masses M > 30A/© for Advanced LIGO, or 
(1 + z)M > 10 5 M© for LISA. As previously reported for 
equal-mass mergers [1, 2] the merger contribution to the 
signal tends to dominate strongly for these larger mass 
systems. For the non-equal-mass cases the merger plays 
a role over similar range of masses, though the level to 
which the merger dominates the overall SNR is signif- 
icantly diminished for very unequal masses (right pan- 
els of Figures 4 and 5) compared to the equal-mass case 
(left, panels; see also the figures presented in [1, 2]). In 
some ways the the equal-mass case is exceptional, rather 
than representative. For observations of IMBH merg- 
ers (M > 100Af©) with Advanced LIGO. however, the 
merger always dominates, as the relatively sharp wall in 
low-frequency sensitivity effectively wipes out the inspi- 
ral contribution. 


IV. WAVEFORM COMPARISON 

Numerical relativity now provides a reasonably clear 
picture of the late stages of merger, specifically in the 
form of the waveforms, as discussed above. While con- 
siderable progress has been made in understanding how 
to detect inspiral signals and characterize how they de- 
pend on system parameters, including mass ratio and 





5 




FIG. 3: Scaling of the Fourier amplitude for different, mass ratios by ■Jfj in the left panel (a.), and by r/ in the right panel (b.), 
which appears to be an excellent approximation. 




FIG. 4: SNR contours for Advanced LIGO with q — 1 (a.) and q = 1/6 (b.). Note, when comparing the two panels, that the 
masses are total masses, which determine the overall waveform amplitude. The solid lines correspond to the SNR calculated 
from the full waveform, including the merger, while the dotted lines correspond to the SNR contribution from the portion of 
the signal with frequency lower than the Schwarzschild ISCO frequency. 


spin, there is little similar work addressing observations 
over the signal-space of mergers. Most observational 
work so far has considered these effectively as unmod- 
eled sources. In [15], we examined the relationships be- 
tween the merger waveforms and the physical motion of 
the source, emphasizing simple common features in order 
to form a general characterization of nonspinning merg- 
ers. These features have observational consequences as 
well. In particular, we noted general similarity in the 
late-time portions of dynamics and waveforms, which we 
now revisit. 


In Fig. 6 we plot model waveforms for the 1:1 case, 
the 4:1 case with its amplitude rescaled with the leading- 
order i] dependence, and the difference 5h = hi — h 2 
between these two waveforms. Bv doing so, we see that 
the apparent phase agreement shown in [15] is, not sur- 
prisingly, partially an artifact of aligning all the phases 


at the peak strain amplitude, and thereby enforcing a 
node in Sh at that time. However, the merger is unique 
in the suppression of the final beat prior to ringdown in 
Sh, and enhancement of that beat if a tr/2 phase shift is 
applied. This can be seen in the time series of Fig. 6, but 
is most evident in the Fourier representation of Fig. 7. 
This extended frequency agreement also provides a sim- 
ple explanation for the rj difference in scaling between 
the inspiral and merger. The amplitude scales with to 
a good approximation for both the inspiral and merger, 
but the time interval spent within a given frequency bin 
scales as tj~ 3 / 8 to leading order for the inspiral, but is 
nearly constant for the merger in the case of moderate 
mass ratios. This results in a relative rj 5 ' 8 amplitude 
scaling of h, and therefore SNR, between the inspiral and 
merger. 

While SNR is certainly the most relevant statistic for 
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FIG. 5: SNR contours for LISA with q — 1 (a.) and q — 1/20 (b.). The solid lines again correspond to the full waveform 
SNR, while the dotted lines correspond to the SNR contribution from frequencies lower than the Schwarzschild ISCO frequency. 
While the observable range of high-SNR mergers is reduced by a factor of several at 20:1 from what was seen in the 1:1 case, 
sources are still easily detectable to large redshift over a similarly broad mass range. 



t(M) 


FIG. 6: Differences in 1:1 and rescaled 4:1 mass ratios. 


detection purposes, it tells us little about the details of 
the waveform, in particular the evolution of the phase, 
which may be critical when answering questions regard- 
ing signal characterization. The “match” [19j is a useful 
statistic for more detailed waveform comparisons, as it is 
sensitive to small differences in waveform phase. For any 
two waveforms, hi and h 2 , the match M is defined using 
the noise-weighted inner product (Eq. (4)): 

M = (7) 

v (hi\hi)(h 2 \h 2 ) 

The match can be viewed as the fraction of the matched- 
filter SNR that is recovered by using h 2 as a filter to 


search for h\, rather than using hi itself (the optimal 
filter). The left panel of Fig. 8 shows a typical com- 
parison for the 6:1 case, which should have the strongest 
higher harmonics among the numerical simulations stud- 
ied here. Also, we show a frequency-based comparison 
for the same ease in the right panel of Fig. 8, where 
the (2, ±2) modes can be seen to dominate the signal 
power until well into ringdown (indicated here by a ver- 
tical dashed line). Nonetheless, it is still possible that a 
sub-dominant mode may modulate the signal to a suf- 
ficient degree to significantly diminish the recoverable 
SNR, or to impact the template member that has the 
highest likelihood, if only the dominant mode is used as 






7 



FIG. 7: The Fourier amplitude of the waveform difference 
in Fig. 6, both with the phase shift shown there, as well as 
shifted by tt/2 to illustrate the dependence of the power at 
high frequencies on the phase alignment. 


a filter. We will investigate this possibility, as well as de- 
velop an appropriate formalism for including all modes 
in a calculation of the match. 

We can further calculate the SNR of the difference in 
waveforms, Sh (see Fig. 9), which is essentially a mea- 
surement of our ability to distinguish two waveforms from 
each other. This simple statistic is related to the “mis- 
match”. 1 — M, as well as the SNR, p (see also [30]): 


(5h\Sh) = (hi - h. 2 \hi - fi 2 ) 

= (h-i\hi) + (h-2\h 2 ) — 2<fti|h a > 


— |ft.2|) + 2\h.i\\h 2 \ ( 1 


(hi\h 2 ) 


hi\\h 2 \ 


( 8 ) 


where \h,i\ = a/ ( fii|fii), and the final approximation 
comes from assuming that the SNRs of hi and h, 2 are ap- 
proximately equal (a similar observation has been made 
in [30]). We note that the curves “1:1”, “rescaled 4:1” 
and “1:1 - 4:1” in Fig. 7, are simply the integrands of iftij, 
h 2 \ and \J (5h\5h), respectively, without noise- weighting. 
The latter is roughly an order of magnitude smaller than 
| fill during the entire merger phase. In this case we can 
see from Eq. (8) that the match among moderate-mass- 
ratio mergers is likely to be quite high for ground-based 
interferometers. We can therefore expect that, for in- 
stance, a small subset of merger waveforms would be ca- 
pable of sufficiently covering a large range of nonspinning 
parameter space for detection purposes, but the apparent 
mass ratio degeneracy in the merger will have a negative 
impact on parameter estimation efforts. 

The SNR has a trivial inverse proportionality with lu- 
minosity distance, so swapping SNR and luminosity dis- 
tance can give you, for instance, the distance horizon at 
which the difference between waveforms can be detected. 


by setting the SNR. at some fixed threshold. We do so in 
Fig. 9, using the same Sh as above as an example. We 
use a fixed SNR of p = 10 as the threshold of detectabil- 
ity. The values in Fig. 9 have an interesting implication 
in light of Eq. (8), in that for sources farther than the 
distance horizon, we cannot distinguish between a 1:1 
waveform at Dl, or a waveform with mass ratio ij at a 
distance (q/0.25 )Dl 


To further investigate the implications of the apparent 
degeneracy in mass ratio for moderate-mass-ratio merg- 
ers, we calculate the match between the full equal-mass 
and 4:1 waveforms, including all available harmonics, as a 
function of the source orientation. In Fig. 10, we show as 
an example a comparison with mass M = 3 x 10 fi M© for 
LISA and M — 100 M© for Advanced LIGO, although the 
result will be qualitatively similar for any masses where 
the merger is emphasized relative to the inspiral in the 
whitened waveform (see the Appendix and Fig. 11) for a 
given detector (i. e. any case where the signal merges at 
or below the peak sensitivity of the detector). Figure 10 
shows a sky map of the matches for Advanced LIGO and 
for LISA, where the match at each point corresponds to 
the equal-mass waveform calculated on that point on the 
source sky. The match is maximized over the orientation 
of the 4:1 waveform in the <3? direction at a given 6, where 
4> is the azimuthal coordinate and © is the polar coor- 
dinate. The maximum at a given © then corresponds to 
the “best” match, and the minimum to the “minimax” 
match 31] . The azimuthal sky position is degenerate 
with the orbital phase, so that the maximization proce- 
dure is identical to finding the maximal value in the az- 
imuthal direction for a particular inclination. We do not 
maximize over the inclination, since the spin-weighted 
spherical harmonics are a more complicated function of 
polar angle. In the sky maps, it is clear that the maxi- 
mization over polar angle would occur at the poles, where 
the quadrupole modes are most dominant. This is consis- 
tent with our previous results showing the striking sim- 
ilarity of the quadrupole radiation across modest mass 
ratios. We observe the expected “north/south” symme- 
try, since all nonspinning waveforms evolve in a fixed 
plane. The azimuthal asymmetry is greatest in the or- 
bital plane, where the fractional luminosity of the higher 
harmonics relative to the dominant quadrupole modes is 
greatest. We note that the sky map woiJd be uniform for 
single mode matches, so the structure in Fig. 10 is the 
result of the harmonic content, and therefore requires the 
formalism contained in the Appendix in order to maxi- 
mize quasi-analytically. The average match over the sky 
of the source for the cases in Fig. 10 is 0.96 for Ad- 
vanced LIGO. and 0.95 for LISA. This means that the 
1:1 waveform can be considered an effective template (in 
the sense of [31]) for typical Advanced LIGO mass ratios 
for a large fraction of source orientations. 




FIG. 8: In the left panel, we show a time series representation of a 6:1 mass-ratio 10 h Mg black-hole binary at a distance of 10 
Gpc, using the waveforms from our numerical simulation. This is a typical example of the (2, +2) modes constituting the vast 
majority of the overall power content of the waveform. This is further demonstrated in the Fourier-series representation (right 
panel), where the (2, +2) modes dominate until well into the ringdown, the onset of which is approximately indicated by the 
dashed vertical line. 




FIG. 9: Luminosity distance where the difference between the 1:1 and rescaled 4:1 mass ratio waveforms is detectable with an 
SNR of 10 for initial LIGO and LISA, also referred to as the distance horizon. This can be interpreted as being the maximum 
distance at which we can distinguish these two sources with each interferometer. 


V. CONCLUSIONS 

We have applied a model for nonspinning late inspiral- 
merger-ringdown waveforms to answer questions regard- 
ing the implications of including the merger phase in data 
analysis efforts. We have verified that, while the merger 
contributes a smaller fraction of the total SNR as we 
deviate from the equal mass case, it still dominates for 
moderate mass ratio, providing nearly the entirety of the 
detectable signal for for ground-based observations of the 
more distantly observable IMBIT systems. In addition, 
we have studied the commonality previously observed in 
the phase evolution of the merger waveform for moder- 
ate mass ratios. While this commonality bodes well for 
detection, since the equal-mass merger waveform alone 


would do well as a search filter for all moderate mass 
ratios, this has negative implications for signal charac- 
terization. Indeed, we have demonstrated that the 1:1 
waveform can be considered an effective template for de- 
tection for a wide range of source orientations for both 
Advanced LTGO and LISA. 
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FIG. 10: Sky map of matches between the equal-mass and 4:1 waveforms for Advanced LIGO (a.) and LISA (b.). Throughout 
the text, <f> is the azimuthal coordinate and © is the polar coordinate. The specific case shown corresponds to a redshifted mass 
of M = 100 M© for Advanced LIGO, and M = 3 x 10 (, M© for LISA (the luminosity distance is irrelevant for this calculation), 
but the result will be qualitatively similar for any mass M > 100M© for Advanced LIGO and M > 2 x 10 B M© for LISA, due 
to the constancy of the detector response and the similar spectral content of the noise over the band of the signal for those 
cases, i. e. cases where the “whitening” procedure (see the Appendix) emphasizes the merger. The sky location corresponds to 
the position on the sky of the equal-mass waveform, with the 4:1 waveform being rotated in the $ direction to maximize the 
match. Therefore, for a fixed 0, the maximum in $ will correspond to the “best” match, and the minimum to the “minimax” 
match. 


Appendix A: Generalized Phase Maximization 

We are interested in generalizing the procedure pre- 
sented in [31] for maximizing the match (Eq. (7)) with 
respect to the initial orbital phase constants between a 
target or exact (label X) waveform and a template or 
approximate (label A) waveform. Specifically, whereas 


the previous method is restricted in its validity to ra- 
diation that is quadrupole-only. we wish to derive the 
general method for maximizing the match for arbitrary 
harmonic content. Wherever possible, we preserve the 
original notation from (31 . 

For the exact and approximate waveforms, we can rep- 
resent the measured strain waveform as 
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h A ' x ( ■ ■■) = F+h + + F x h x = !R [Fe ,K h] 
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where Af m x = coa(m(p A X ), \f; x = sin (m<p A - x ), h. Um = 
Fjbr ra |cos(ro</>), form = F|h« m | sin(m^), Fe tK = F + + 
?F X is the complex beam pattern function, and 0 and o 
describe the angular position on the source’s sky (with 
k absorbed into the definition of <j>). Since <j> can be ab- 
sorbed into <p, the following procedure maximizes over 
the relative azimuthal orientation as well as the orbital 
phase. For this analysis, we assume a common source 
polar angle 9 for the exact and approximate waveforms, 
although the procedure could be further generalized to al- 
low maximization/ minimization over all relevant angles. 

As in [31], we wish to find the “best/' and “minimax” 
match. We therefore wish to form an appropriate basis 
in which we can decompose the exact and approximate 
waveforms separately, and subsequently find the projec- 
tion of the resulting approximate “vector” on the exact 
“vector”. Conceptually, in [31] the procedure amounted 
to finding the ellipse resulting from projecting the circle 
that the approximate waveform makes in its 2-plane onto 
the 2-plane formed from the decomposition of the exact 
waveform, where the 2-planes are the spaces spanned by 
the orthonormal bases constructed using the exact and 
the approximate waveforms. In our analysis, we extend 
this concept to folding the minimum and maximum ra- 
dius resulting from finding the sum of projections of ap- 
proximate circles for the available modes on the exact 
planes corresponding to those same modes. We could 
alternatively include in the sum the cross contributions 
from particular approximate circles for a given mode on 
the planes of all available exact modes. However, while 
such a result would be more directly related to the SNR 
achievable by using the approximate waveform as a tem- 
plate, the inclusion of cross-mode contributions would 
be unphysical and less useful as a gauge for potential pa- 
rameter estimation. We therefore include only like-mode 
contributions, although the following derivation can be 
trivially altered to include all cross-mode contributions, 
and the final result will be the same in all but the most 
exotic cases. 

To form the desired bases, we first construct a pair of 
“whitened” vectors [32], as shown in Fig. 11, in both 
the approximate and exact planes, to account for the 


detector response to the raw 

2), 


h A ' X ' 

"'iiZm 



cc 


-%2tx ft 
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where “h” denotes the Fourier transform of h. With these 
whitened vectors, the noise-weighted inner product (4) 
can be easily calculated in the time domain: 

pOG 

(h\ ho) = / dth!*{t)h! 2 {t). (A3) 

J o 

Instead of attempting to construct a single orthonor- 
mal basis, we generate an orthogonal (not normal) basis 
for each im mode, with the normalization chosen so that 
the sum over modes is normalized, i.e. 


A,X 

e iZm 
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(A4) 


This expression yields an appropriate normalization over 
all modes, since e AX = YLim e nim is normalized by con- 
struction, with the individual im modes being appropri- 
ately weighted by their relative barycentric power and 
the response of the detector. (A4) is therefore a set of 
orthogonal bases which are all constrained by the total 
signal power, and by their common dependence on the 
orbital phases of the exact and approximate waveforms. 
We therefore retain the original two degrees of freedom 
as in [31]. 

If we focus only on like-mode contributions to 
the match, we can construct the projection operator, 
Px( e a ) i all d the resulting projection p a , of a vector e A 
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FIG. 11: Examples of “whitened" waveforms [32] that we use to form a basis for calculating the match. Examples for Advanced 
LIGO are shown in the left panel, and examples for LISA are shown in the right panel. The ordinate values are arbitrarily 
scaled. 


onto the X-plane, 

Pc = Px (e£) = ]T X (m > (A5) 
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defined as 

e iim = cos(ma) ef e * + sin (ma) . (A6) 

where a is an arbitrary initial orbital angle a for the 
approximate waveform. Substituting for from (A6) 
into (A5) yields 

Pc = 22 [cos(ma) p Um + sin(ma) p 2tm \ , (A7) 
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While (A9) is trivial to maximize or minimize analyt- 
ically for the case of a single mode as in [31], the case of 
multiple modes generally requires a numerical solution. 
However, if we assume a single mode (or mode pair) is 
significantly larger than any other mode, then we can 
specify an approximate solution for the value of a that 
yields the ‘‘best’' match. Generally, the condition for ex- 
tremizing (A9) is given by 


v , / 

Pntm = Pxif-nlrn) — / , ~ 


,A 

'n£m 


[zt (4mi e L) 


(A8) 


y: [m ( Atm - Bfm) sin(2ma) 

t.m 

— 2m,Cim cos(2ma)j =0. (All) 


If we again focus only on like- mode contributions for sim- 
plicity, then |p Q [ can be expressed in a form which rep- 
resents, geometrically, a sum of ellipses, given by 

\p a \ 2 = [Atm cos 2 (mo) + Bem. sin 2 (ma) 

f,rn 

+'2Ct;rn cos(ma) sin(ma)| , (A9) 


We can then apply the aforementioned assumption that 
a single mode pair dominates. In geometric terms, this 
means that we assume that the semi-major axis for the 
dominant mode(s) in Eq. (A9) is larger than the quadra- 
ture sum of the semi-major axes of all other modes. In 
this case, the largest value for Eq. (A9) will occur very 
near the a that maximizes the dominant mode(s), with 
the other modes providing at most a small perturbation. 



12 


This condition can be expressed as 
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where £,V1 corresponds to the dominant mode, and we 
only include the larger roots of Eq. (A9) corresponding 
to the “best” match. For all cases in this paper, £ = 
\M\ = 2, with IPcCmI = \PaC(.-M)\ by symmetry, so 
that the condition for (A12) in this case will be tm =/= 
£|A4|. Finally, we can calculate the condition on a for 
maximizing the match under these assumptions: 


C^best 


2 M 



(Acm - Bcm) 
cm — Bcm) 2 + 4C£ 



We reiterate that Eq. (A13) is not valid if subdominant 
modes contain comparable power to the dominant mode 
or mode pair, and a similar method cannot be used to 
find the “minimax” a. In these cases, Eq. (A9) can only 
be solved numerically. Even if a numerical solution is 
required, this method is still more efficient than a brute- 
force maximization over <p A and ip x , as it makes it a one- 
dimensional search over a. Since A( m , Be m ■ and Ct m are 
all less than unity, the error in the match will be of the 
same order as the sampling interval in a over the range 
(0, 2tt], assuming modes with very large m are negligible. 
For instance, in this work we only include t<4 modes, 


with m = + 4 the largest relevant m mode. Therefore, in 
order to calculate the match to three significant digits, 
we take 10 4 5 6 7 8 9 samples of a and record the global max- 
imum and minimum, corresponding to the “best” and 
“minimax” matches, respectively. One could implement 
a more clever algorithm, such as Brent’s method, if the 
required accuracy for the match makes the sampling pro- 
cedure too computationally expensive. We have verified 
that the minimum and maximum from Eq. (A9) agrees 
with the minimum and maximum found using a Nelder- 
Mead simplex over the two-dimensional c p A -tp x space. 

We note that, even if we were to include the cross- mode 
contributions in our derivation, the result would remain 
a sum of a set, albeit a much larger set, of ellipses. It 
should be noted that the exclusion of cross-mode con- 
tent makes this representation of generalized matches less 
closely related to a matched-filtered SNR calculation or 
even, potentially, to the maximum likelihood estimator, 
in that it does not account for circumstances where the 
maximum likelihood occurs at an incorrect value for the 
parameters. In such cases, this method will instead ig- 
nore all maxima except the local maximum determined 
from the largest match of like-modes. In most cases, the 
contribution of cross-mode terms to any match calcula- 
tion will be negligible compared to like-mode contribu- 
tions, so this local maximum will also be the global maxi- 
mum. In that case, the optimized parameter choice found 
from this procedure will be consistent with the maximum 
likelihood value, and the resulting match will be a true 
representation of the fraction of recoverable SNR from 
matched filtering. 
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